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Abstract 

In this paper, we are concerned with global existence and optimal decay rates of solutions for 
the three-dimensional compressible Hall-MHD equations. First, we prove the global existence 
of strong solutions by the standard energy method under the condition that the initial data are 
close to the constant equilibrium state in iJ 2 -framework. Second, optimal decay rates of strong 
solutions in L-norm are obtained if the initial data belong to L 1 additionally. Finally, we 
apply Fourier splitting method by Schonbek [Arch.Rational Mech. Anal. 88 (1985)] to establish 
optimal decay rates for higher order spatial derivatives of classical solutions in IJ-framework, 
which improves the work of Fan et al. [Nonlinear Anal. Real World Appl. 22 (2015)]. 
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1 Introduction 


The application of Hall-magnetohydrodynamics system (in short, Hall-MHD) covers a very wide 
range of physical objects, for example, magnetic reconnection in space plasmas, star formulation, 
neutron stars, and geodynamo, refer to |l|-|6j] and the references therein. Recently, Acheritogaray et 
al. [3] derived the Hall-MHD equations from the two-fluid Euler-Maxwell system for electrons and 
ions through a set of scaling limits or from the kinetic equations by taking macroscopic quantities in 
the equations under some closure assumptions. They also established the global existence of weak 
solutions for periodic boundary condition. In this paper, we investigate the following compressible 
Hall-MHD equations in three dimensional whole space R 3 (see 0]): 


' p t + di v(pu) = 0, 

(pu)t + di v(pu <g> u) — pAu — (p + i/)Vdmt + 'VP(p) = (curLB) x B, 

(curlR) x B 


B t — curl(u x B) + curl 


( 1 . 1 ) 


= A B, divR = 0, 


where the functions p , u, and B represent density, velocity, and magnetic field respectively. The 
pressure P(p) is a smooth function in a neighborhood of 1 with P'{ 1) = 1. The constants p and u 
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denote the viscosity coefficients of the flow and satisfy physical condition as follows 


p > 0, 2p + 3zz > 0. 

To complete the system m, the initial data are given by 

(p,u,B)(x,t )| t=0 = (p 0 (x),u 0 (x), B 0 (x)). 
Furthermore, as the space variable tends to infinity, we assume 

lim (p 0 - 1 ,u 0 ,B 0 )(x) = 0. 

|x|—>-oo 


( 1 . 2 ) 


(1.3) 


Obviously, the compressible Hall-MHD equations transform into the well-known compressible MHD 
equations when the Hall effect term curl ^ (curiB) x B ^ j s neglected. 

When the density is constant, Chae et al.[8] proved local existence of smooth solutions for large 
data and global smooth solutions for small data in three dimensional whole space. They also showed 
a Liouville theorem for the stationary solutions. Chae and Lee jfl] established an optimal blow-up 
criterion for classical solutions and proved two global-in-time existence results of classical solutions 
for small initial data, the smallness conditions of which are given by the suitable Sobolev and Besov 
norms respectively. Later, Fan et al.[lo| also established some new regularity criteria, which also 
are built for density-dependent incompressible Hall-MHD equations with positive initial density by 
Fan and Ozawa 11]. Maicon and Lucas [12J] proved a stability theorem for global large solutions 
under a suitable integrable hypothesis and constructed a special large solution by assuming the 
condition of curl-free magnetic fields. Fan et al. [13| established the global well-posedness of the 
axisymmetric solutions. Chae and Schonbek [1J] established temporal decay estimates for weak 


solutions and obtained algebraic time decay for higher order Sobolev norms of small initial data 
solutions as follows 


I V k u 


3+2 fc 


L 2 + ||V fc H(t)|| L 2 <C(l + t)-— , he N 


for all t > T*(T* is a positive constant). Furthermore, Weng [15(] extended this result by providing 
upper and lower bounds on the decay of higher order derivatives. For the compressible Hall-MHD 
equations D, Fan et al.[16| proved the local existence of strong solutions with positive initial 
density and global small solutions (classical solutions) with small initial perturbation. They also 
established optimal time decay rate for classical solutions as follows 


\\(p-l,u,B)(t )\\ L 2 < C(l + t) 4. 


(1.4) 


Here, they required the initial perturbation is small in Lf 3 -norm and bounded in L 1 -norm. 

Recently, the study of decay rates for solutions to the MHD equations has aroused many re¬ 
searchers’ interest. First of all, under the lL 3 -framework, Li and Yu [l7j] and Chen and Tan 18( not 
only established the global existence of classical solutions, but also obtained the time decay rates 
for the three-dimensional compressible MHD equations by assuming the initial data belong to L 1 
and L q (q € [l, |)) respectively. More precisely, Chen and Tan 18f] built the time decay rates 


\v k (p-i,u,B)(t)\\ H3 - k <c(i + ty 


(1.5) 


where k = 0,1. The time decay rates has also been established by Li and Yu jl7[ for the case 
q = 1. Motivated by the work of Guo and Wang [191 ], Tan and Wang [20[] established the optimal 
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time decay rates for the higher order spatial derivatives of solutions if the initial perturbation 
belongs to H v n H~ s (N >3 ,sG [0, §)). More precisely, they built the following time decay rates 


l|V fe (p - l,u,B)(t)\\ HN - k < C(l+t)-^, 


where k = 0,1,..., N — 1. Motivated by the work 22], we (see 23]) establish the following time 
decay rates for all t > T*(T* is a positive constant), 


\7 k (p — l)(t)\\ H 3-k + WV ku (t)\\ H 3-k < C(1 + t) 4 : 
\7 m B(t)\\ H 3-m < C( 1 + t)-^, 


where k = 0,1, 2, and m = 0,1,2, 3. It is easy to see that the time decay rates (11.61) is better than 
decay rates m since (11.61) provides faster time decay rates for the higher order spatial derivatives 
of solutions. 

In this paper, we hope to establish the global existence and time decay rates of solutions for 
the compressible Hall-MHD equations (ll.ll) - (ll.3l) . First of all, we construct the global existence of 
strong solutions by the standard energy method under the condition that the initial data are close 
to the constant equilibrium state (1,0,0) in H 2 - norm. Second, if the initial data in L 1 — norm are 
finite additionally, the optimal time decay rates of strong solutions are established by the method 
of Green function. Precisely, we obtain the following time decay rates for all t > 0, 

3+2 k 

\\(p — l)(t)\\ H 2-k + \\u(t)\\ H 2-k + \\B(t)\\ H 2-k < C(1 + t) 4 , 


where k = 0,1. This framework of time convergence rates for compressible flows has been applied to 
other compressible models, refer to [24-1271]. Although magnetic held equations (II.IIP are nonlinear 
parabolic equations, we hope to establish optimal time decay rates for the second order spatial 
derivatives of magnetic held under the condition of small initial perturbation. In order to achieve 
this goal, we move the nonlinear terms to the right hand side of (11. IIP and deal with the nonlinear 
terms as external force with the property on fast time decay rates. Then, the application of Fourier 
splitting method by Schonbek [2lJ helps us to establish optimal time decay rate for the second order 
spatial derivatives of magnetic held as follows 


|V 2 H(t)|| i2 <C(l + t)-4. 


Finally, one focus on establishing optimal time decay rates for higher order spatial derivatives of 
classical solutions to compressible Hall-MHD equations. More precisely, we prove that the global 
classical solution ( p,u,B ) of Cauchy problem (ll.ll) - (ll.3l) has the time decay rates (11.61) . Obviously, 
these time decay rates improve the results (11.41) by Fan et al.[16| since we build faster time decay 
rates for higher order spatial derivatives of classical solutions. 

Notation: In this paper, we use H s (M?)(s € R) to denote the usual Sobolev spaces with norm 
|| • || h 3 and L p (R 3 )(1 < p < oo) to denote the usual L p spaces with norm || • ||^p. The symbol \7 l 
with an integer l > 0 stands for the usual any spatial derivatives of order l. For example, we dehne 


V k v = {d%Vi\ |o| = k, i = 1,2,3}, v = (ui,u 2 ,u 3 ). 


We also denote ^(f) ■= f ■ The notation a <b means that a < Cb for a universal constant C > 0 
independent of time t. The notation a ~ b means a < b and b < a. For the sake of simplicity, we 
write f fdx := J r3 fdx. 
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First of all, we establish the global existence and optimal decay rates of strong solutions for the 
compressible Hall-MHD equations (ll.ll) - (|1.3l) . 

Theorem 1.1. Assume that the initial data (po — 1, no, -Bo) £ R 2 and there exists a small constant 
do > 0 such that 

IK A) - 1,uo,B 0 )\\ H 2 < So, 

then the problem (HU-H3D admits a unique global strong solution (p, n, B) satisfying for all t > 0, 

rt 


\\(P- hu,B)(t)f H2 + [ (||Vp(s)||^ 1 + ||V(n, B){s)\\ 2 H2 )ds < C ||( Po - l,«o, 
Jo 


B, 


0 / 11 ' 


Furthermore, if ||(po — 1,no, Bo)||i,i is finite additionally, then the global strong solution ( p,u,B) 
has following decay rates for all t > 0, 

l|V fc (P - l)(f)||jf 2 -k + \\V k u(t)\\ H 2 -k < (7(1 + t) 4 , ^ 

\\V m B(t)\\ H 2-m ^Cil + t)-^ 1 , 
where k = 0,1, and m = 0,1, 2. 

Remark 1.1. For any 2 < p < 6, by virtue of Theorem U. 1\ and the Sobolev interpolation inequality, 
we obtain time decay rates as follows 

,-4(1- 


II (p ~ 1)(*)I|lp + ll u (*)l|i> < + t) 2 ' 

||V fc B(t)|| LP < C(1 + 

where k = 0,1. Furthermore, in the same manner, we also have 

|| (/i — l)(^)l|i/°° + ||n(t)||i,°° < C(1 + 1) 4 , 

\\B(t)\\ L °° <C(l + t)~l 

Second, we build time decay rates for the time derivatives of global strong solutions. 

Theorem 1.2. Under all the assumptions in Theorem II. 11 the global strong solution ( p,u,B ) of 
Cauchy problem (fTTP-tfOD has the decay rates 

\\pt(t)\\ H i + ||nt(i )|| i 2 < C( 1 + t) 4 , 

\\B t {t)\\ L 2 < C(l + t)~l 

for all t > 0. 

Furthermore, we establish optimal decay rates for the higher order spatial derivatives of classical 
solutions to the compressible Hall-MHD equations. 

Theorem 1.3. Assume that the initial data (po — 1,uq,Bq) £ H 3 n L 1 and there exists a small 
constant Eq > 0 such that 

II (Po - 1,^0iB 0 )|| H 3 < £q, (1.8) 

then the global classical solution ( p,u,B) of the problem fjl.ljl has the time decay rates 

WV k (p — l)(t)\\ H 3-k + \\V k u(t)\\ H 3-k < C(1 +1) 4 , 

||V m B(t)||^3-m < C( 1 + t 

where k = 0,1,2, and m = 0,1,2, 3. 


(1.9) 
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Remark 1.2. Compared with the decay rates of linearized systems of (ED stated in Proposition 
ED EHD gives optimal decay rates of the solutions and its spatial derivatives (except for the third 
order spatial derivatives of density and velocity) in L 2 -norm to the nonlinear problem (ll.l|) ~ (ll.3h . 
Here the decay rate of solutions to nonlinear system is optimal in the sense that it coincides with 
the rate of solutions to the linearized systems. 

Remark 1.3. By virtue of the Sobolev inequality and the results EHD in Theorem El then the 
global classical solution ( p,u,B ) has the time decay rates 

\\(p — l)(t)\\LP + \\u(t)\\LP < C(1 + t) 

||v fe RCO|| LP + 

where k = 0,1, and p € [2, oo]. Hence, the decay rate of classical solution ( p,u,B ) converging to 

3 

the equilibrium state (1,0,0) in L°°-norm is (1 + t)~ 2 . 

Remark 1.4. It is easy to see that EHD provides faster time decay rates for higher order spatial 
derivatives of global classical solutions than El- Hence, the results in Theorem El improve the 
work of Fan et al. & 

Remark 1.5. Although we only established the time decay rates under the -framework in The¬ 
orem El the method here can be applied to the H n (N > 3) -framework just following the idea as 
Gao et alS]. Hence, if (po — 1,uq,Bq) € H N n L 1 (A^ > 3), then the global solution (p,u,B) has 
the time decay rates 

l|V fe (p — l)(£)||#jv-fc + ||V fc 'u(t)||^jv-fc < C(1 + t) 4 , 

||V m R(t)||^JV-m < C( 1 + t 

where k = 0,1, ...,N — 1, and m = 0,1,2, 

Finally, we build decay rates for the mixed space-time derivatives of global classical solutions. 

Theorem 1.4. Under all the assumptions in Theorem El the global classical solution (p,u,B) of 
the problem El-El satisfies the time decay rates 

\\V k p t (t)\\ H2 - k + l|V fe ^(t)|| L 2 < C(l + t)^, 

\\V k B t (t)\\ L2 <C(l + t)- 7 -^, 

where k = 0,1. 

This paper is organized as follows. In section 2, we establish some energy estimates that will 
play an essential role for us to construct the global existence of strong solutions. Then, we close the 
estimates by the standard continuity argument and the global existence of strong solutions follows 
immediately. Furthermore, we build the time decay rates by taking the method of Green function 
and establish optimal time decay rates for the second order spatial derivatives of magnetic field. 
Finally, we also study decay rates for the time derivatives of density, velocity and magnetic field. 
In section 3, we establish the optimal decay rates for the higher order spatial derivatives of global 
classical solutions and mixed space-time derivatives of solutions. 
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2 Proof of Theorem 11.11 and Theorem 11.2 


In this section, we will establish global existence and optimal time decay rates of strong solutions 
for the compressible Hall-MHD equations. Indeed, computing directly, it is easy to deduce 

1. 


(curlS) x S = (S • V)S - ^V(|S| 2 ), 


and 


curl (a x B) = u(divS) — (u ■ V)S + (S ■ V)u — S(divtt). 
Then, denoting g = p — 1, we rewrite (11.11) in the perturbation form as 

' g t + divu = Si, 

< ut — pAu — (p + z/)Vdmz + = S 2 , 

. B t — AB = S 3 , divS = 0, 

where the function S*(i = 1, 2, 3) is defined as 

" Si = — gdivu — u -Vg, 

S 2 = —u • V u —h(g)[pAu + (p + z^Vdivu] — f(g)Vg + g(g) 


B ■ VS - -V(|S| 2 ) 


S 3 = —u ■ VS + S • Vu — Bdivu — curl 
Here the nonlinear function of g is defined by 


9{q) ( S • VS — - V(|S| 2 ) 


h{g) = 


f(Q) = P ' [6 + l) - 1 , g(g) = 1 


f? + 1 g 1 ^+1 

The initial data are given as 

(g,u,B)(x,t)\ t=0 = (g 0 , u 0 , B 0 )(x) -A- (0,0,0) as |x| -> 00 . 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 


2.1. Energy estimates 

First of all, suppose there exists a small positive constant 5 satisfying following estimate 

\\{g,u,B)(t )\\ H 2 := || 0 (t)||tf 2 + |K *)|| ff 2 + \\B(t )\\ H 2 < 6, (2.5) 

which, together with Sobolev inequality, yields directly 

1 3 
- < g+ 1 < -. 

2 “2 

Hence, we immediately have 

\f(g)\,\h(g)\ < C\g\ and {g^^(g)\,\h (k) (g)\,\f (k) (g)\ < C for any k>l, ( 2 . 6 ) 

which will be used frequently to derive a priori estimates. 

We state the classical Sobolev interpolation of the Gagliardo-Nirenberg inequality, refer to 
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Lemma 2.1. Let 0 < m, a < l and the function f G then we have 

I|V q /||lp < ||V™/||^||V'/||£ 2 , 

where 0 < 0 < 1 and a satisfy 

First of all, we will derive following energy estimates. 

Lemma 2.2. Under the condition (USD, then for k = 0, 1, we have 

j t \\V k (g,u,B)\\l 2 +C\\V k+ \u,B)\\l 2 < 5\\\7 k+1 g\\ 2 L2 . 


(2.7) 


( 2 . 8 ) 


Proof. Taking k- th spatial derivatives to (I2.1l) i and m 2 respectively, multiplying the resulting 
identities by S7 k g and V k u respectively and integrating over M 3 (by parts), it is easy to obtain 


~ /d V ^| 2 +\^ k u\ 2 )dx + J{ft\S7 k+l u\ 2 + (ji + ^)|V fc divu| 2 )dx 


/ 


= / V k S i • V k g dx + / V k S 2 • V k u dx. 


S 


(2.9) 


Taking k- th spatial derivatives to (| 2 . 1 |) ^, multiplying the resulting identity by \7 k B and integrating 
over M 3 (by parts), we have 


1 d_ 

2 dt 


\V k B\ 2 dx + J \X7 k+1 B\ 2 dx = j V k S 3 -V k Bdx. 

Adding (12.91) to (12.101) . it follows immediately 

~J(\V k g\ 2 +\V k u\ 2 + \S/ kB \ 2 ) dx +f{n\V k+1 u\ 2 + (R+u)\W k dw u \ 2 +\V k+1 B\ 2 )da 
= / V fc 5i • V k g dx + [ V k S 2 ■ V k u dx + I V k S 3 ■ V k B dx. 


( 2 . 10 ) 


( 2 . 11 ) 


1 d_ 

2 dt 


For the case k = 0, then the differential identity (12.111) has the following form 

f (\g\ 2 +\u\ 2 + \B\ 2 )dx + j (fi\Vu\ 2 + (// + v)\divu\ 2 + \VB\ 2 )dx 
= j Si ■ g dx + j S 2 ■ u dx + j S 3 ■ B dx = I\ + I 2 + I 3 . 

Applying (12.5p . Holder, Sobolev and Young inequalities, it is easy to obtain 

h < II^IIl 3 l|divu||^211^11^6 + HelMlVelMMUe 
< llellHi||v«|| La ||V0|| L 2 + ||e||Hi||Ve|| La ||v«|| La 


( 2 . 12 ) 


(2.13) 


5(\\\7g\\ 2 L2 + \\Vu\\ 2 L2 ). 
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Integrating by parts and applying (12.61) . Holder, Sobolev and Young inequalities, it arrives at directly 

— J h(g)(fiAu + {fj, + v)V&\vu)udx 
~ / (h'(g)Vg ■ u + h{g)Vu)S7udx 


< 


< 


< 


llVell^llltlUellVtilliS + llelUacUVullia 
(llell^ + lIVtillHOdlVellia + IIVullia) 
^llVellia + UVullia). 


(2.14) 


Hence, with the help of (12. 6 p . (|2. 14[) . Holder, Sobolev and Young inequalities, we deduce 


h < (\\u\\ L s\\Vu \\ L 2 + ||e||L3||Ve|| L 2 + \\g(g)\\ L ^\\B\\ L3 \\VB\\ L 2 )\\u\\ L e 

+ <S(l|Vf?||la + ||Vu||£a) 

< (||«[|hi||V«|| l2 + Hell/rillVelUa + ||S||„i||VS|| i2 )||Vu|| L2 (2.15) 

+ <y(l|Ve||ia + ||V«||ia) 

< ^(||V^||| 2 + ||V«||ia + ||VS|||a). 

Integrating by part and applying (12.61) . Holder and Sobolev inequalities, it arrives at 


— J curl \g(g)(B ■ VS)] Bdx 

= ~ j g(e){B ■ VS)curLBete (2.1.6) 

<||^)|| L c»||S|| L oc||VS|| L 2 ||curlS || L2 

^ \\ B \\hA\^ B \\ 2 L 2. 

Hence, with the help of (12.161) . Holder, Sobolev and Young inequalities, we deduce 

h < (N| L 3||VS|| i2 + ||S|| L 3||V«|| i a)||S|| i 6 + ||5|| H a||VB||£a 

< (\\u\\ H i + l|5|| H i)(||V«||ia + ||VS||| 2 ) + <5||VS||| 2 (2.17) 

< ^(HV< 2 + ||VS||ia). 


Substituting (12.131) . (12.151) and (12.171) into (12.121) and applying the smallness of 5, it arrives at 
directly 

j t f(\ Q ?+\u\ 2 + \B\ 2 )dx + j(g\Vu\ 2 + \VB\ 2 )dx<6\\Vg\\ 2 L2 . (2.18) 

For the case k = 1, then the differential identity (12.111) has the following form 


~ /(l V ^| 2 +l Vw l 2 + |VS| 2 )dx + J(g\V 2 u\ 2 + (// + i/)|Vdivu | 2 + |V 2 S| 2 )d a 
= I VSi-Vgdx+ [ V5 2 • \7u dx + [ VS 3 - VS dx = Ih + II 2 + il 3 . 


(2.19) 


Applying Holder, Sobolev and Young inequalities, we obtain 


n i < (II^IIl 3 11divit11x,® + ||n|| L 3||V0|| i 6)||V 2 ^|| L2 
<(ll^ + IHkO(llv 2 5||l 2 + ||v 2 < 2 ) 
<- 5 (l|v 2 e|l ! 2 + l!v 2 < 2 ). 
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Similarly, it is easy to deduce 


Ih < (|M| L 3||V«|| L « + ||^)||l-||V 2 u|| l2 )||V 2 u|| l2 

+ (ll/(e)llL»l|Ve|| L 6 + lb(0)llL-||s||L3||VB|| L .)||v 2 «|| L2 
< (llellfl, + \\u\\ H i + \\B\\ Hl )(\\V 2 0 \\l 2 + ||V 2 U||| 2 + ||V 2 H|| 2 2 ) 
<<5(l|V 2 ^||i 2 + ||V 2 u||i 2 + ||V 2 i?||| 2 ). 


Integrating by part and applying (12.61) . Holder and Sobolev inequalities, we obtain 


- j Vcur \[g(g)(B • VH)]VH dx 
= - J V[g(g)(B • VH)]VcurLB dx 

< (||V 5 (e)IU8||B||i8||VS|U« + ^( e )|| L =c||VB|| L 3||VS|| L «)||VcurlS|| L2 
+ ||»(e)|U« \\B\\l~ ||V 2 J B|| i2 ||VcurlH|| i2 

< (||V 2 e || L2 ||VH|| i2 + ||V5|| i3 + ||H|| h2 )||V 2 H|| 2 2 . 


( 2 . 22 ) 


Applying (12.221) . Holder, Sobolev and Young inequalities, it arrives at directly 


Ih < (||u|| L 3||VH|| L 6 + || J B|| L 3|iVu|| i 6)||V 2 H|| L2 

+ (|!V 2 £|| i2 ||VH|| L2 + ||VH|| L 3 + ||H|| h2 )||V 2 H|| 2 2 (2.23) 

< <5(||V 2 ti|| 2 2 + ||V 2 H|| 2 2 ). 


Substituting (12.201) . (12.21 j) and (12.231) into (12.191) . then we obtain 

j t |(|V £ »| 2 +|Vu| 2 + |V J B| 2 )dx + y‘(^|V 2 u| 2 + |V 2 H| 2 )dx<«5||V 2 ^||i 2 , 
which, together with (12.181) . completes the proof of the lemma. □ 


Next, we derive the second type of energy estimates involving the higher-order spatial derivatives 
of g and u. 


Lemma 2.3. Under the condition (|2.5I) . then we have 


d 


j t \\V\g lUl B)\\l 2 +C\\V 3 (u,B)\\ 2 L2 <d\\V 2 e\\ 2 L 2. 


(2.24) 


Proof. Taking k = 2 specially in (12. lip , we deduce immediately 

^/(| v2 ^| 2 + l v 2 u | 2 + | v 2 H| 2 )dx+/(^| v 3 u | 2 + (^+i/)| v 2 d ivu | 2 + | v 3 H| 2 )dx 
= [ V 2 Si • \7 2 q dx + [ V 2 S 2 ■ V 2 u dx + [ V 2 5 3 • V 2 H dx. 


(2.25) 
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Applying Holder, Sobolev and Young inequalities, it is easy to obtain 
— J SJ 2 {gd\vu)V 2 gdx 

= — J (V 2 £divu + 2V^Vdivu + gSJ 2 d\vu)\7 2 gdx 


< 


(\\Vu\\ L ~\\\7 2 g\\ L 2 + \\Vg\\ L3 \\X7 2 u\\ L6 


\ l ~\\V 3 u\\ L 2)\\\7 2 q\\ L 2 


< (\\v 2 u\\l 2 \\v 3 u\\Uv 2 e\\L* + ||v^|| H i||v 3 n|| L2 + \\Q\\ H 4v 3 u\\Lz)\\v 2 e\\L* 


< 


V 2 u\\l 2 \\v 2 0 \\l 2 + UVellHi + \\q\\h*W 2 q\\1. + ||v 3 u|| 2 2 


<5(||v 2 £ »||i 2 + ||v 3 «||i 2 ). 

Integrating by part and applying Holder, Sobolev and Young inequalities, it arrives at 
- [\7 2 (u ■ Vg)V 2 g dx 

1 „ „ 1 

dx 

< (||v 2 «|U«||v e |U3 + ||Vu|| Lt x J ||v 2 e || L2 )||v 2 f? || L 2 


— (V 2 uV£ + 2V'uV 2 £>)V 2 £ + ^|V 2 f?| 2 divrt 


< 


l|V^|ki||V 2 ^|| L 2||V 3 u|| L 2 + ||V 2 u||£ 2 ||V 3 «||£ 2 ||V 2 £>||£ 2 
< (||V 2 u||| 2 ||V 2 e||| 2 + ||V^|| //1 )(||V 2 e||i2 + ||V 3 u||| 2 ) 
<-5(I|v 2 ^ii!2 + iiv 3 w || 2 2 ). 


The combination of (|2.26|) and (12.271) gives rise to 

J V 2 Si • V 2 ^ dx < 5(\\\7 2 q\\ 2 l2 + ||V 3 u||i 2 ). 


(2.26) 


(2.27) 


(2.28) 


Now, we give the estimate for the second term on the right hand side of (12.251) . By virtue of Holder 
and Sobolev inequalities, we have 


J V(u • \7u)X7Audx 
= J (VuVu + u\7 2 u)\7Audx 


< l|V'u|| L 3||Vu|| L 6||V 3 u|| i 2 + |M| L 3||V 2 u|| L 6||V 3 u|| L 2 

< 11^111,2 l|V 3 w|| ^21| V^ll^a II V 3 xt||2 2 II V 3 xt|| Z/2 + IMI/p II V 3 u||| 2 

<^I|v 3 u|| 2 2. 

In view of (12.61) . Holder and Sobolev inequalities, we have 


(2.29) 


J V(/i(g)(/zAit + (n + v)\7 divu))X7 Audx 

< (l|V/ i (e)|U3||V 2 n|| i6 + \\h(g)\\ L ~\\V 3 u\\ L2 )\\V 3 u \\ L 2 (2 . 30) 

< (l|V^|kl||V 3 W || L 2 + || e || H 2||V 3 «|| L 2)||V 3 u|| i 2 
<^I|v 3 u||| 2 , 


to 
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J V(f(g)Vg)\/Audx 

<(I|V^|||4 + ||/(^)||l-||V 2 0|| l 2 )||V\|| l2 

< (||vf ^|| L2 ||V 2 e || L2 + ||^|| //2 ||V 2 e || i2 )||V\|| i2 

< (l|V^|| L2 + ||e|b 2 )||V 2 f? || i2 ||V 3 u|| L2 

<-5(iiv 2 ^iii 2 + nv 3 n|ii 2 )- 


Similarly, it is easy to deduce 


V 


g(o){B • V-B — V(-|-B| 2 )) 


VA udx 


< (l|V^)|| L 6||S|| L 8||VS|U6 + II^IIl-IIVSII^IIVSH^JIIV^II^ 
+ \\ 9 ( q )\\ l ~\\b\\ L 3\\v 2 b\\ l 4v 3 u\\ L 2 


< (||V 0 || i 6||V J B|| i2 ||V 2 i?|| L2 + || J B||| 2 ||V 3 i?||| 2 ||Vi?||| 2 ||V 3 J B||| 2 )||V 3 u|| i2 

+ \\B\\ h1 \\V 3 B\\ l2 \\V 3 u\\l* 

<<5(||V 2 e||| 2 + ||V 3 u||| 2 + ||V 3 J B||| 2 ). 


(2.32) 


By virtue of the estimates (|2.29l) - (|2.32l) . we obtain immediately 

J V 2 5 2 • V 2 u dx < <5(||V 2 f?|| 2 2 + ||V 3 u|| 2 2 + ||V 3 H|| 2 2 ). (2.33) 

Integrating by part and applying Holder and Sobolev inequalities, it arrives at 


j V 2 (m-VB)V 2 B dx 

< (\\Vu\\ L 3\\\7B\\ L e + |M| L 3||V 2 H|| L6 )||V 3 £|| L2 

< (ll^ll| 2 ||V 3 u||| 2 ||VH||| 2 ||V 3 i?||| 2 + ||u|| h i||V 3 j B|| l 2 )||V 3 H|| l2 

< <5(||V 3 ^||i 2 + ||V 3 J B||| 2 ). 


Similarly, it is easy to obtain 


J V 2 (H • Vu + Hdivu)V 2 £ dx 

< (||V«||i3||VS|| L 8 + || J B|| i 3||V 2 u|| i 6)||V 3 J B|| L2 

< (IKIli 2 ||V 3 u||| 2 ||VH||| 2 ||V 3 H||| 2 + || j B|Ui||V 3 u|| l 2 )||V 3 H|| l2 
<d(||V 3 < 2 + ||V 3 J B||| 2 ). 


ll 
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Integrating by part and applying (12.61) . Holder and Sobolev inequalities, we obtain 


/ 


— / V^curl 


g(g)(B-VB-V(-\B\ 2 ] 


= - / V" 


g(g)(B-VB-V(-\B\ 


V 2 B dx 
V 2 cur LB dx 


< 


\V 2 g(g)\\ L 2 \\B\\ LO o IIVSIUco + ||V 5 (f?)|| L 6 ||VB|| l6 ||VH|| i6 )||V 2 curlB|| i2 


+ (I|v^)IIl6|| j b||l6||v 2 j b|| L 6 + 


L~\\VB\\ L3 \\V 2 B\\ L e)\\V 2 cmlB \\ L 2 


+ \\g(g)\\ L oo\\B\\ L ^\\V 3 B\\ L2 \\W 2 cuv\B\\^ 

< (IHV^I 2 + \V 2 g\\\ L 2 \\B\\ L ~\\VB\\ L ~ + ||V e || L6 ||VH|| L6 ||V J B|| L6 )||V 3 J B|| i2 
+ (l|Ve|| L 6||H|| L6 ||V 2 J B|| L 6 + ||VH|| L 3||V 2 H|| i6 + ||H||^||V 3 J B|| L2 )||V 3 B|| i2 


< (l|V^||| 2 ||V 2 ^||| 2 + ||V 2 ^|| l2 )||B|^ 2 ||V 2 H||| 2 ||V 3 j B||| 
+ ||VB||^ 1 ||V 2 p|| l2 ||V 3 j B|| l2 + \\VB\\ h i\\V 3 B\\ 2 l2 
<-5(l|V 2 ^||i 2 + ||V 3 H||i 2 ). 


(2.36) 


In view of the estimates (12.341) - (12.361) . we obtain directly 


J V 2 B 3 • V 2 B dx < <5(||V 2 £|| 2 2 + ||V 3 n||| 2 + ||V 3 ^||| 2 )- 


(2.37) 


Substituting (I2.28p . (I2.33P and (12.371) into (I2.25D . then we have 

j t J(\V 2 g\ 2 +\\7 2 u\ 2 + \V 2 B\ 2 )dx + J(g\V 3 u\ 2 + \V 3 B\ 2 )dx<5\\X7 2 g\\ 2 L2 , 

which completes the proof of the lemma. □ 


Finally, we will use the equations (12.11) to recover the dissipation estimate for g. 

Lemma 2.4. Under the condition ()2.5I) . then for k = 0,1. we have 

J X7 k u • V k+1 gdx + C\\V k+1 g\\ 2 L2 < || V fc+1 u|| 2 2 + \\V k+2 u\\ 2 L2 + ||V fc+2 H|| 2 2 . (2.38) 

Proof. Taking fc-th spatial derivatives to the second equation of (12.11) . multiplying by V k+1 g and 
integrating over M 3 , then we obtain 


j \7 k u t -\7 k+1 gdx + J \V k+1 g\ 2 dx 

= J S7 k [gAu + (ji + i/)Vdi vu]\7 k+1 gdx + j \7 k S 2 ■ V k+1 gdx. 


(2.39) 


In order to deal with f V k UfV k+1 gdx, we turn the time derivatives of velocity to the density. Then, 
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; equation (12.11b . 

V k u t 

■ \7 k+1 gdx 

If 

V k u • V fc+1 

dt J 



\7 k u • V fc+1 

dt J 


±f 

V k u • V fc+1 

dt J 


±f 

V k u • V fc+1 

dt J 



7&+1, 


(2.40) 


(2.41) 


(2.42) 


Substituting (I2.40P into (I2.39p . it is easy to deduce 

j t J V k u-V k+1 gdx + J \V k+1 g\ 2 dx 

= j \V k divu\ 2 dx + J V^’divu • V fc div(uu)^' + J V k S 2 ■ V fc+1 gdx 
+J \7 k [gAu + (n + v) V divu] V fc+1 gdx. 

For the case k = 0, then applying Holder, Sobolev and Young inequalities, we obtain 

J divu • di v(gu)dx < ||£»||l°° || Vu||| 2 + \\u\\ L 3 ||divu|| L 6 ||Vp|| L2 

~ (ll^lln 2 + IMIfrOGI^IIl 2 + l|V 2 ^||| 2 ) 

<S(\\Vg \\ 2 L2 + \\V 2 u\\l 2 ). 

By virtue of (I2.6p and Holder inequality, it is easy to deduce 

J S 2 ■ Vgdx < (||tx|| L 3||V«|| L 6 + ||^||L-||V 2 u|| L2 )||V^|| i2 

+ (||elU«»||Ve|| L2 + ||5(e)||L-||5|!L3||VH|| L 6)||v^|| L2 
<<5(l|Ve||i 2 + ||V 2 u||i 2 + ||V 2 H||| 2 ), 

and 

/ [gAu + (g + i/)Vdivu]V gdx < ||V 2 u||| 2 + 4Ve\\b- 

The combination of (J2.42I) . (1 2.43 p and (1 2.44 p helps us complete the proof to (|2.38p for the case of 
k = 0. As for the case k = 1, applying Holder, Sobolev and Young inequalities, we deduce 

[ Vdivu • \7dlv(gu)dx 

(2-45) 


(2.43) 


(2.44) 


^ (l|V^|| L 3||divu|| L 6 + ||g|| Z/ oc||Vdivu|| L2 )||V 2 u|| L2 
+ (||V£|| L 3 ||Vu || l6 + |MM|v 2 e|| L 2 )||v 2 u || i2 
<5(l|V 2 ^||i 2 + ||V 2 u||| 2 ). 


With the help of Holder inequality and Lemma 12.31 it arrives at 

j VS 2 ■ V 2 gdx < 5(\\V 2 g\\ 2 L2 + ||V 3 ^||| 2 + ||V 3 S||£ 2 ), 


(2.46) 
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and 

/ V[//A u + (/x + i/)Vdivu]V 2 gda; < ||V 3 w||| / 2 + e||V 2 e|| 2 2 . (2.47) 

The combination of (12.451) . (j2.46l) and ()2.47l) gives rise to the proof of (|2.38l) for the case of k = 1. □ 

2.2. Global existence of strong solutions 

In this subsection, we shall combine the energy estimates that we have derived in the previous 
section to prove the global existence of strong solutions. Summing up (12.81) from k = l (l = 0,1) to 
k = 1, then we obtain 


d_ 

dt 


+C\\V(Vu,VB)r Hl _i <<511^+^11^ 


which, together with (I2.24p . then it arrives at 
d 


dt 


V\g,u,B)\\ 2 H2 _ l +C\\V l+ \u,B)\\ 2 H2 _ l < < 5C 1 ||V ,+1 0|& 1 _ I . 


(2.48) 


On the other hand, summing ()2.38l) from k = l (l = 0,1) to k = 1, we obtain immediately 

j Y / • v fc+ Vx + c 2 ||v ,+ 1 e|& 1 _ l < c (liv'+^ii^-, + ||v^ 2 5 ||^_ ( ) . (2.49) 

at l<k<I'' 

Multiplying (|2.49l) by 25C\/C2 and adding the resulting inequality to (|2.481) . then it arrives at 

; £?(t) + C 3 (||V ,+1 <>|&i-« + \\V l+1 (u,B)\\ 2 H2 ^ < 0. (2.50) 

where £f ( t ) is defined as 


d_ 

dt 


-2/-J.N _ DM|2 1 25Cl 

C 2 


£f(t) = \\V(g,u, B)\\ 2 h2 _i + 


Y [ V k u-V k+1 gdx. 

s'US'•' 


Kk< 1 


(2.51) 


By virtue of the smallness of 5, it is easy to obtain 

Ol V Z G?, U, B)\\ 2 H2 _ t < £ 2 (t) < C4V l (g, u, B) ||^ 2 _*. 

Choosing l = 0 in (12.501) . integrating over [0, t] and applying the equivalent relation (12.511) . we obtain 

\\(g,u, B)(t )\\ H 2 < C\\(g 0 ,u 0 ,B 0 )\\ H2 . 


Then, by the standard continuity argument (see Theorem 7.1 on page 100 in [29j]), we close the 
estimate (12.51) . Thus, we extend the local strong solutions to be global one and the uniqueness of 
global strong solutions is guaranteed by the uniqueness of local solutions that has been prove by Fan 
et al. [l6j. Therefore, choosing l = 0 in (12.50(1 . integrating over [0,t] and applying the equivalent 
relation (12.5ip . we obtain 

\\(g,u,B)(t)\\ 2 H2 + [\\\Vg(T)\\ 2 Hl + \\(Vu,VB)(T)\\ 2 H2 )dT < C\\(q 0 ,u 0 ,B 0 )\\ 2 h2 , (2.52) 

Jo 

which completes the proof of the global existence of strong solutions. 
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2.3. Decay rates of strong solution 


In this section, we will establish optimal decay rates for the compressible Hall-MHD equations 
(ll.ll) - (ll.3l) . If the initial perturbation belongs to L 1 additionally, we apply the method of Green 
function to establish optimal time decay rates for the global strong solutions. Furthermore, the 
application of Fourier splitting method by Schonbek 2l|] helps us to build optimal time decay rates 
for the second order spatial derivatives of magnetic field. 

First of all, let us to consider the following linearized systems 


' Q t + divw = 0, 

< ut — rAu — (ft + z/)Vdivu + = 0, (2.53) 

. Bt — AB = 0, 


with the initial data 


(g,u,B)(x,t)\ t=0 = (q 0 ,u 0 ,B 0 )(x) -j- (0,0,0) as \x\ oo. (2.54) 

Obviously, the solution ( g,u,B ) of the linearized problem (12.53l) - (12.54l) can be expressed as 


(g,u,B) tr = G(t ) * ( g 0 ,u 0 ,B 0 ) tr ,t > 0. 


(2.55) 


Here G(t) := G(x,t ) is the Green matrix for the systems (12.531) and the exact expression of the 
Fourier transform G(£, t ) of Green function G(x, t ) as 


G(£,t) = 


\+e x - t -X-e x + t 

A_|_—A_ 

—i£(e x + t —e x — t ) A_|_e A + f —A_e 


-i^t( e x +t_ e x- t) 


A+-A- 

0 


A+-A- 


A_i 

HP 


A+-A_ 


D Aot 


hx3 — 


w 


0 

0 

e Alt / 3x 3 


where 


Ao = ~r\£\ 2 , Ai = -|£| 2 , 


X + = -U + -u)\f\ 2 + iJ\^- U+-u) |e| 


A- = - ICI 2 - »W l£l 2 - + 2^J Kl 4 - 

Since the systems (12.571) is an independent coupling of the classical linearized Navier-Stokes equa¬ 
tions and heat equation, the representation of Green function G(£, t ) is easy to verify. Furthermore, 
we have the following decay rates for the systems (I2.53I) - (I2.54I) . refer to [30]. 


Proposition 2.5. Assume that ( g,u,B ) is the solution of the linearized systems (j2.53p ~ (|2.54p with 
the initial data (go,uo, Bo) € L 1 D H 2 , then 

||V^||| 2 <G(||(^o^o)|||i + ||V fe (^o^o)||| 2 )(l + rK 
\\V k uf L2 < C (||(eo,«o)|lii + \\V k (go,uo)\\ 2 L2 ) (1 + t)~l~ k , 

||V fc H||2 2 < C (\\Bo\\ 2 l1 + ||V fc Ho||i 2 ) (1 + t)~l~ k 

for 0 < k < 2. 
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In the sequel, we want to verify some estimates that play an important role for us to derive 
decay rates for the compressible Hall-MHD equations f)2. 1 li - (I2.4ji . 

11(51,52,53)11x1 < 5{ \\S7q\\ L 2 + ||V«|| H 1 + HVSlI^i), 

||(5 1 ,5 2 ,5 3 )||x 2 < <5(||Ve|| L2 + \\Vu\\hi + llVH^i), (2.56) 

||V(5 1 ,5 2 ,5 3 )||x2 < S(\\V 2 q\\l* + l|V 2 u|| i2 + \\S7 2 B\\ L 2) + ||V(e,B)|| H i||V 2 (u,5)|| H i. 

Now, we establish the decay rates for the compressible Hall-MHD equations (|2.ip ~ (|2.4|) . 

Lemma 2.6. Under the assumptions of Theorem \l.ll the global strong solution ( g,u,B) of problem 
(EH)-(Ha]) has the time decay rates 

\\V k g(t)f H2 _ k + \\'V k u(t)\\ 2 H 2 -k + \\\7 k B(t)\\ 2 H 2 -k < C(1 + t)~i~ k (2.57) 


for k = 0,1. 


Proof. Taking l = 1 specially in (|2.50l) . it arrives at directly 

j t P 2 (t) + C 3 (||V 2 £»Hi 2 + ||V 2 u||^ + \\V 2 B\\ 2 m ) < 0, (2.58) 

where S 2 (t) is defined as 

£i(t) = \\Vg\\ 2 m + \\Vu\\ 2 Hl + \\VB\\ 2 h1 + ^jvu- V 2 gdx. 

With the help of Young inequality and smallness of 5 , it is easy to deduce 

C^\\V(g,u,B)\\ 2 Hl < Slit) < C 4 \\V(g,u,B)\\ 2 Hl . (2.59) 

Adding both hand sides of (12.581) by ||V(^, u, B)\\ 2 l2 and applying the equivalent relation (|2.59[) . 
then we have 

j t S 2 (t) + CS 2 (t)<\\V(g,u,B)\\ 2 L2 . 

In view of the Gronwall inequality, it follows immediately 

S 2 it) < SH0)e~ ct + f e- c ^\\S7ig,u,B)iT)\\ 2 L2 dT. (2.60) 

Jo 

In order to derive the time decay rate for S 2 (t), we need to control the term || V(g, u, B)\\ 2 l2 . In 
fact, by Duhamel principle, we can represent the solutions for the problem (12.1jl - ()2.4l) as 

(f?, U, B) tr it) = Git) * (go, uo, B 0 ) tr + f G(t - s) * (5i, 5 2 , S 3 ) tr (s)ds. (2.61) 

Jo 


Denoting 

E(t)= sup (l + r)S(||Ve(r)||| rl + ||V«(r)||i rl + ||V5(r)||| rl ) > 

0 <T<t 
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which, together with (12.61 p . (|2.56D and Proposition 12.51 gives directly 

/i i 

IIV(e, U, B)(t) \\ L 2 < C(l + t)-t + C / (H(5i, 5 2 , 5 3 )(t)|| l i + || V(Si, S 2 , S 3 )(t)|| L 2) (1 + t - t)~Ut 

Jo 

pi 

< C{\ + +c 6 (II V 0 (r)|| H i + ||Vu(r)|| H i + ||V5(t)|| h i) (1 + t - t)~Ut 

Jo 

+ C f II V(e, B)(t)\\ H 1 || v 2 (+ 5)(t)|| h i (1 + 1 - t)-Ut 
J o 

< C7(l + t)“t + . 


(1 + t — t) 4(1 + 7 -) *dr 


0 


+ Cy/E(i) 


_ 1 


(1+t — r) 2(1 + r) 2 dr 
5 


l‘ \\\/ 2 (u,B)(r)\\ 2 m dr 
.Jo 


< C( 1 + 1 )~4 + C<5/E(t)(l + i) -3 

<(l+t)-i(l + <5 v / £(*)), 


where we have used the fact 

rt 


[ (l + t-rni+r)' 

./o 


r dr 


+ / (l+t-T)- r (l+T)- r (iT 


< 1 + 



(1 + r) r dr + ( 1 + 


—r r t 



(1 + t — r) T ’dr 

/ jo \ j\ 

<(l + t)~ r , 

for r = | and r = | respectively. Thus, we have the estimate 

IIV(e,«, fi)(i)Ilia < C( 1 + t)-i(1 + <5£(t)). 

Inserting (J2.62I! into (12.601! . it follows immediately 

£f(t) <£f(0)e~ ct + C [ e- c{t ~ T) (l + T)- 5 2 (l + 5E(T))dT 

Jo 

<£f(0)e~ ct + C(l + 6E(t)) [ e -C(t-r)(i + T )-ldr 

Jo 

< £ 2 (0)e~ Ct + C( 1 + <SE(t))(l + t)“5 
<C(l + ^(t))(l + t)-f, 
where we have used the fact 

rt 


f e- C(t - T) (l+r)-tdr 

Jo 

= So + J* e_C(t ~ r)(1 + r )"^ T 


< e ^ 


P(1 + r)' 

J 0 


2 dr + 1 + 


-C(t-r ) d7 


(2.62) 


(2.63) 


< C(l + t)“5 . 
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Hence, by virtue of the definition of E(t) and (I2.63|) . it follows immediately 

E(t) <C(l + 5E(t)), 
which, in view of the smallness of <5, gives 

E(t) < C. 


Therefore, we have the following time decay rates 

l|V(?(£)llfl- 1 + ||Vu(i)||k + \\VB(t)\\ 2 H2 < C(l + t)~l. (2.64) 

On the other hand, by (12.611) . (12.561) . (|2.64l) and Proposition 12.51 it is easy to deduce 

\\(g,u,B)\\ 2 L2 <C(l + t)--2+C [\\\(S 1 ,S 2 ,S 3 )\\l 1 + \\(S 1 ,S 2 ,S 3 )\\ 2 L2 )(l + t-r)-ldT 

Jo 

< C( 1 + 1 )~§ + C f 5 (|| V<?(t )||£ 2 + ||Vu(r)||^ + || V^t)^) (1 +1 - t)~Ut 

Jo 

< 0(1 + *)“§ +C [ (l+t-r)“i(l + T) _ idr 

Jo 

< 0 ( 1 + *)“§, 


where we have used the fact 

(1 + t- r) _ t(l + r) _ idr < 0(1 + t)~T 

Hence, we have the following decay rate 

IkWIll 2 + II^WIIl 2 + ll-^WIli 2 ^ C(1 + 0 2 - (2.65) 

Therefore, the combination of (12.641) and (12.651) completes the proof of the lemma. □ 



Finally, we establish optimal decay rates for the second order derivatives of magnetic field. 


Lemma 2.7. 

decay rate 


Under the assumptions of Theorem 11.11 then the magnetic field has the following 


V 2 B(t)\\ L 2<C(l + t)~l 


( 2 . 66 ) 


Proof. Taking k = 2 in (|2.10l) . it follows immediately 

\V 2 B\ 2 dx + J \V 3 B\ 2 dx = J V 2 5 3 • V 2 H dx. 

By Holder, Sobolev and Young inequalities, we obtain 


(2.67) 


J V 2 (H • Vu + Bdivu)\7 2 B dx 

< (||Vu|| L 3||VH|| L 6 + ||H|| L ^||V 2 u|| L2 )||V 3 H|| i2 (2 . 68) 

< (||V«|M|V 2 H|| l2 + ||VH|| H1 ||V 2 u|| i2 )||V 3 H|| L2 
<||V(u,H)||^||V 2 (u,H)||i 2 +5||V 3 H||| 2 . 


18 

















Global Existence and Optimal Decay Rates of Solutions for Compressible Hall-MHD Equations 


It follows from (|2.34l) and (12.361) that 

J V 2 (u ■ VB)\7 2 B dx < \\Vu\\ 2 h1 \\V 2 B\\ 2 l2 + 6\\V 3 B\\ 2 l2 , 


(2.69) 


and 


g(e) [b ■ VB - V(^|-B| 2 ) 


V 2 B dx 


— j V curl 

<\\W 2 0\\UWB\\ 2 Hl +6\\W 3 B\\l 2 . 

Substituting (|2.68p ~ (12.70l) into (12.671) and applying the time decay rates (12.571) . then we obtain 

j t J \V 2 B\ 2 dx + J \V 3 B\ 2 dx 

< \\W(u,B)\\ 2 Hl \\V 2 (u,B)\\l 2 + \\V 2 q\\ 2 l 2 \\VB\\ 2 h1 

< (1 + t)~i{ 1 + t)~I + (1 + t)~ 1(1 + i)-i 

<(i+r 5 - 


(2.70) 


For some constant R defined below, denoting the time sphere(see 21]) 

/■ i \ 

R 




1 + t 


it follows immediately 


\V 3 B\ 2 dx> [ |£| 6 |I3| 2 d£ 

Jr 3 /So 


> 


[ I«I 4 |S| 2 <« 

1 + t JR 3 /S 0 


> 


R 

1 + t 


/ lel 

J R 3 


4l B\ 2 d£ — 


R 

1 + t 


'So 


lei Wdt, 


or equivalently 


\V 3 B\ 2 dx > I \V 2 B\ 2 dx - ( 

1 + t JR 3 V 


R 

1 + t 


WBrdx. 


The combination of (12.711) . (12.721) and (12.571) yields directly 

4 [ \V 2 B\ 2 dx H-—— f |V 2 S' 2 

dtj 1 1 l + tj 1 


dx 


< 


16 


a +ty 


\VB\ 2 dx + C(1 + t) 


-5 


< 


(l+t)~ 2 (l+t)-2 + (l+t) 


-5 


< C(1 + f) 2 , 

where we have chosen R = 4 in (12.721) . Multiplying (12.731) by (1 + f) 4 , we obtain 

±{(l + t) 4 \\V 2 B\\ 2 L2 } <C(l + t)-l 


(2.71) 


(2.72) 


(2.73) 


(2.74) 
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Integrating (12.7411 over [0 ,t], then we have the following decay rate 

\\V 2 B(t)\\ 2 L2 <C(l + t)-l. 

Therefore, we complete the proof of the lemma. □ 

Proof of Theorem ll.lt With the help of (12.5211 . Lemma l2.6l and Lemma [2.71 we complete the 
proof of Theorem 11.11 


\Bt\\v = 


A B — u • VB + B ■ \7u — Bdivu — curl 


g( 0 )[B.X7B-±X7(\B\ 2 ) 


L 2 


~ l|AB|| L 2 + ||«||l 3||VB|| L 6 + ||B|| L 3||Vn|| L6 + ||Vg(^)|| L 2||B|| L 6||VB|| L 6 
+ || ff (2)|| L oo||VB|| i3 ||VB|| L 6 + Me)\\ L ~\\B\\ L ~\\V 2 B\\ L 2 

< ||V 2 B|| L 2 + \\(u,B)\\ m \\V 2 (u,B )\\ L 2 + ||V^|| L 2 ||VB||^ + ||VB||^i||V 2 B|| L 2 

< {l + t)~i + (1+ f)-i(l+f)"t + (l + t)-l(l + t)-l 

< C(l + t)~v. 

In view of the decay rates (I2.76|l - (l2.79p . we complete the proof of the lemma. 


(2.75) 


(2.76) 


2-4- Proof of Theorem \1.2\ 

In this subsection, we establish the decay rates for the time derivatives of strong solutions. 

Lemma 2.8. Under the assumptions of Theorem M.il the global strong solution ( g,u,B ) of problem 
(GHi-dMD satisfies 

11 I Iff 1 + IKWIIl 2 < C( 1 + t) 

\\B t {t)\\ L 2 < C(l + t) 4 . 

Proof. By virtue of the equation ()2.1H i and decay rates (HD, we have 

Weth 2 = ||divu + £»div« + u ■ X7g \\ L 2 

< ||divn|| i2 + ||^|| LO o||divu || L 2 + ||Vp|| i3 \\u\\ L 6 

< c{i + t)~u 

Similarly, it follows immediately 

||^7 gt 11/^2 = ||Vdivtz + V^divu + ^Vdivn + Vn - Vg + u ■ V 2 ^||^2 

< ||Vdivu || L 2 + ||V^|| L 3||V'u|| L 6 + ||(£»,u)|| L oo||V 2 (£>,u )|| L 2 

< ||V 2 u|| L 2 + ||(£, u)\\ H 2 ||V 2 (£», u)\\ L 2 

< C(l + t)~i. 

By virtue of the equation (12. ID ?, decay rates (11.711 and estimate (I2.56p . we have 

IMIl 2 = 11/uA u + in + z^Vdfvu - Vg + S 2 \\l 2 

< ||V 2 n|| L2 + ||Ve || L 2 +<J(||Ve || L 2 + ||V(u,B)|| h i) 

<(l + t)-f + (l + t)-f 

< C(1 +t) _ l 

By virtue of (I2.ip 2. (11.71) . Holder and Sobolev inequalities, we obtain 


(2.77) 


(2.78) 


(2.79) 


□ 
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Proof of Theorem 11.21 With the help of Lemma 12.81 we complete the proof of Theorem 11.21 


3 Proof of Theorem 11.31 and Theorem 11.41 

In this section, we first establish optimal time decay rates for the higher order spatial derivatives 
of global classical solutions under the condition of small initial perturbation in iL 3 -norm and finite 
initial perturbation in L'-norm. Furthermore, we also study the decay rates for the mixed space- 
time derivatives of global classical solutions. 

First of all, Fan et al.(see (3.2) on Page 430 in JL0]) have established following estimate 

\\(g,u,B)(t)\\ H 3 < C\\(qo,u 0 ,B 0 )\\ H 3 < Ce 0 . (3.1) 

Thus, the inequality (12.6p also holds on under the condition of (11.81) . 


3.1. Proof of Theorem 11.31 

Just following the idea as Lemma 12.61 it is easy to establish optimal decay rates for the global 
classical solutions. For the sake of brevity, we only state the results in the following lemma. 

Lemma 3.1. Under the assumptions of Theorem 1 1 . 31 the global classical solution ( g,u,B ) of prob¬ 
lem (I2.1D - (I2.4I) satisfies for all t > 0, 

l|V*e(t)|&3- h + ||V fe u(f)||^3- fc + ||V fe 5(l)||^-* < C(l + t)-l- k , (3.2) 


where k = 0,1. 


Next, we establish optimal time decay rates for the second order spatial derivatives of magnetic 
field and enhance the time decay rates for the third order spatial derivatives of magnetic field. 

Lemma 3.2. Under the assumptions of Theorem 1 1. 31 then the magnetic field has following decay 
rate for all t > 0, 


||V 2 5(1)||^ <C(l + l)-4. 


(3.3) 


Proof. Taking k = 3 in (12.101) . it follows immediately 
j \V 3 B\ 2 dx + j |V 4 5| 2 dx 


1 d 

2 dt 


= / V 3 


— u • VB + B ■ V« — Bdivu — curl ( g(g) ( B ■ WB — -V(|5| 2 ) 


V s Bdx 


(3.4) 


= IIIi + III 2 + III 3 + Hh- 


By virtue of (13.11) . Holder and Sobolev inequalities, it arrives at 

IIIi = j (V 2 wV5 + 2VuV 2 B + u\7 3 B)\7 A Bdx 

< (HV^MlVHlUe + ||Vu|| L 3 ||V 2 5|| L 6 + ||n|| i 3||V 3 5|| L 6)||V 4 5|| i2 

< ||V 2 u|| 2 3 ||V 2 B|| 2 2 + ||Vw|| 2 3 ||V 3 5|| 2 2 + (e + e 0 )\\V 4 B\\ 2 l2 

< ||Vu||^ 2 ||V 2 5||^ 1 + ( e + e o)||V 4 5|| 2 2 . 
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In view of the Sobolev and Young inequalities, we obtain 

III 2 = ~ J(V 2 B\7u + 2 VBV 2 m + B'V 3 u)'V 4 Bdx 

< (HV^IMlVullis + ||VB|| L 6||V 2 u|| L 3 + ||-B||l°°||V 3 u|| L 2)||V 4 -B||i2 (3 . 6) 

< (||V 3 i?|| i 21|Vii||jji + ||V 2 S|| L2 ||V 2 U || H 1 + ||V J B|| Hl ||V 3 «|| L2 )||V 4 J B|| i2 

< ||V 2 B|| 2 , 1 ||V«|| 2 , 2 + IIVSII^IIV^II 2 , + e[|V 4 5|| 2 a . 

In the same manner, we get 

Ilh < l|V 2 5||^||Vu||^ + ||VS|| 2 rl ||V 3 «|| 2 a + e||V 4 5|| 2 a . (3.7) 

Applying (12.6p . (13.11) . Holder, Sobolev and Young inequalities, it is easy to deduce 

III 4 < (\\V 3 g(g)\\ L 4B\\ L ~\\VB\\ L ~ + ||V 2 ^)|| i 3||V J B|| i6 ||VH|| iC x>)||V 4 J B|| i2 
+ (l|V 2 ^)llL3p||L-||V 2 H|| L 6 + ||V^)|| L 6||VH|| i 6||V 2 H|| L6 )||V 4 H|| L2 
+ (II^ 7 5 , (^)IIl 6 II^IIl 6 I|V 3 H|| L 6 + ||fl , (^)IU oo l|VH|| Z/ 3||V 3 H|| L 6)||V 4 H|| i2 
+ \\ 9 (q)\\l~\\b\\ l ~\\v 4 b\\1 2 

< ||VB||^ \\V 2 B\\ 2 Hl + ||V 2 A?||| 2 IIV 2 H||^ + II VB\\ 2 h1 IIV 3 H|| 2 2 + (e + e 0 )IIV 4 A?||| 2 

< ||VH|^ 1 ||V 2 J B||| 1 + ( £ + eo )||V 4 i?||| 2 . 

Substituting (13.5D - (13.8D into (13.41) and applying the smallness of s and £q, it is easy to deduce 

j t j\V 3 B\ 2 dx + j |V 4 J B| 2 dx<||V(u,H)||| 2 ||V 2 (rr, J B)||^ 1 , (3.9) 

which, together with the time decay rates (13.21) . yields directly 

j t J\V 3 B\ 2 dx + J \X7 4 B\ 2 dx < (1 + t)~ 5 . (3.10) 

Similar to (12.721) . it is easy to deduce 

J\V 4 B\ 2 dx >Y^~ t J\V 3 B\ 2 dx - j \V 2 B\ 2 dx. (3.11) 

The combination of (13.21) . (13.101) and (13.111) gives directly 

3/ |V 3 B|2<;i + _1_/ |V 3 B|2(;i 

~ (1 -I- i)2 J |V 2 -B| 2 ete + (1 + t) 5 

<(l + t)- 2 (l + t)-f+ (l + t)“ 5 

<(1 + t)S, 

which, together with (12.731) . yields directly 

J t j (|V 2 5| 2 + |V 3 H| 2 )dx + Y^- t J (|V 2 H| 2 + |V 3 H| 2 )dx < (1 + t)-l (3.12) 
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Multiplying (|3.12p by (1 + t) 4 , it arrives at 

— [(1 + t) 4 ||V 2 i?||^-i] < (7(1 + t) 2, 
which, integrating over [0,f], gives 

||V 2 S(t)||2 fl <C(l + t)~l. 

Therefore, we complete the proof of the lemma. □ 


In order to establish optimal decay rate for the third order spatial derivatives of magnetic field, 
we need to improve the decay rate for the second and third order spatial derivatives of velocity. 
Indeed, following the idea as the compressible MHD equations(see |23|), it is easy to verify the time 
decay rates (13.211 also holds on for k = 2. For the convenience of readers, we also introduce the 
method to improve the decay rates for the second order spatial derivatives of density and velocity 
here. 


Lemma 3.3. Under the assumptions of Theorem \1.3l the global classical solution (g, u, B ) of Cauchy 
problem (12.111 - (12.41) has 


-||V 2 (£, u)f Hl + H|v 3 u||^i < c 5 [(1 + 1 )- 5 + s 0 \\v 3 e\\ 2 L2 ]. 

Proof. Taking k = 2 specially in (12.91) . then we obtain 

~ /(l v2 e| 2 +l v2 u| 2 )dx + J(g\V 3 u\ 2 + (g + v)\\7 2 dlvu\ 2 )dx 
= [ V 2 Si • V 2 g dx+ [ V 2 S 2 ■ V 2 U dx. 


(3.13) 


(3.14) 


Integrating by part and applying (13.21) . Holder, Sobolev and Young inequalities, we obtain 


/ V 2 Si • V 2 p dx = j V(gdivu + u ■Vg) ■ V ' 3 gdx 

< (||Ve|| £ 3||V«|| L 8 + ||V 2 u|| L 3||^|| L 6 + \\V 2 g\\ L3 \\u\\ L e)\\\/ 3 g \\ L 2 

< (l + t)- 5 + e||V 3 e||| 2 . 


From the estimates (12.2911 and (|2.30l) . it is obtain the estimates 

J V 2 [— u-Vu — h(g)[p,Au + (g + u)V d\vu}\\/ 2 udx < £o||V 3 tt|| 2 2 . (3.16) 

Integrating by part and applying (|2.6I) . (1 3.2 [) . Holder and Sobolev inequalities, we obtain 

f V 2 [-f(g)\7g]\7 2 u dx 

< (ll/(e)IU-l|v 2 ^|| L 2 + ||v/( e )|| i 3||v^|| L6 )||v 3 u|| L2 

< (Heb- + ||v f? || i 3)||v 2 e || i2 ||v 3 u|| i2 (3 ‘ 17) 

<l|v^||| 1 ||v 2 ^||i 2 + e ||v 3 u||i 2 

<(l + t)- 5 + e||V 3 u||| 2 . 
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In the same manner, it is easy to deduce 


1. 


V 2 g(g) ( B ■ WB — -V(|i?| 2 ) 


V 2 u dx 

(Vg{g)BVB + g{g)VBVB + g{g)B\7 2 B)V 3 udx 


< (I|V£'||l 6 || j B|| L 6 ||VS || L 6 + 

+ \\9(q)\\l^ \\B\\ L 6\\\ /2 B\\t3\\'V 3 u\ 


L°° ||VB|| i 3||V5|| L 6)||V 3 u|| L 2 


(3.18) 


Il3||v u\\ L 2 

< ||V 2 0 || 2 2 ||V1?|| 2 2 ||V 2 j B|| 2 2 + ||VS||| 3 ||V 2 S||| 2 
+ ||VH|| 2 2 ||V 2 H|| 2 3 +£||V 3 u|| 2 2 
<(l + t)- 5 + e||V 3 u||i 2 . 

In view of the estimates (13.161) — (13.181) . it is easy deduce 


/ 


V"S 2 • V"u dx < (1 + 1 )~ 5 + eoll'V 3 u||i 2 


7 3„.||2 


(3.19) 


Substituting (13.151) and (13.191) into (13.141) and applying the smallness of e and eoi we obtain 
d 


dt 


Taking k = 3 in 
1 d 


{\S7 2 g\ 2 + \S/ 2 u\ 2 )dx + ii / jV 3 u| 2 dx< (l + t)- 5 + e||V 3 f ? ||i 2 


specially, then we have 


(3.20) 


~ y(|V 3 ,| 2 + |V 3 u| 2 )dx + J (g\V 4 u\ 2 + (g + v)\\7 3 divu\ 2 )dx 
= j V 3 (— ^divu — u • Vf?)V 3 gdx + J V 3 [— u- Vu — h(g)(fj,Au + (/i + i/)Vdivtt)] V 3 u dx 


+ / V d 


-/(f?)Ve + g(e)(B- VB --V(|H| 2 )) 


(3.21) 


V 3 u dx 


= /Vi + IV 2 + IV 3 + IV 4 + + /^6 + /Vr- 


Applying the Holer and Sobolev inequalities, it is easy to deduce 

IV-y < (||V 3 e|| i2 ||Vu|| L ~ + ||V 2 ^|| i 6||V 2 u|| L 3)||V 3 ^|| L2 

+ (||V^|| L 3||V 3 u|| i6 + || e ||^||V 4 u|| L2 )||V 3 ^|| L2 (3.22) 

<eo(||V 3 ^|| 2 2 + ||V 4 u|| 2 2 ). 

Similarly, it is easy to deduce 

IV 2 <e 0 (\\V 3 g\\ 2 L2 + \\V 4 u\\ 2 L2 ). (3.23) 

Integrating by part and applying decay rates (13.2 f) . Holder and Sobolev inequalities, it arrives at 

IV 3 = j V 2 (u • Vu)V 4 u dx 

< (||Vu|| L 3||V 2 u|| L 6 + ||u|| i 3||V 3 W || L 6)||V 4 u|U 2 (3.24) 

<||Vu|| 2 3||V\||i 2 + ( £ + eo )||V 4 u||| 2 

<(l + t) 5 + (e + e 0 )||V 4 t(||| 2 . 
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In view of (12.61) . (J3.2I) . Holer and Sobolev inequalities, we obtain 

IV A « J V 2 (h(g)V 2 u)V*u dx 

= J {V 2 h{g)V 2 u + 2 \7h(g)V 3 u + h(g)V A u)V A udx 

< (||V^||i 6 ||V 2 u|| L6 + ||V 2 e|| L 3||V 2 u|| L6 )||V 4 u|| L2 (3.25) 

+ (||V^|| L 3||V 3 u|| L 6 + ||Me)IM|V 4 u|| i2 )||V 4 u|| L2 

< ||V 2 ^|| 2 fl ||V 3 u|| 2 2 + ( e + eo)||V 4 u|| 2 2 
<(l + t) 5 + (e + e 0 )||V 4 tt|| 2 2 . 

Similarly, it is easy to deduce immediately 

IV ,5 = J ( f(g)V 3 g + 2 Vf(g)V 2 g + V 2 / (*?) V g)V A udx 

< (II/(i?)IIl-||V 3 0|| l2 + \\Vq\\ 2 l 6 \\Vq \\ l6 + ||V 2 e || L 3||V^|| i6 )||V 4 u|| i2 (3 . 26) 

< l|V^||^ 1 ||V 3 ^||| 2 + ||V 2 ,||| 2 + \\V 2 g\\UV 2 g\\ 2 L2 +e\\S/ 4 u\\h 

< (1 +1) 5 + e||V 4 u||^ 2 . 

Integrating by part and applying (12.61) . (13.2[) . Holder and Sobolev inequalities, we get 

IVto = ~ J {V 2 g(g)B\7B + 2 \7g(g)V(B\7B) + g(g)V 2 {BVB))V 4 udx 

< (l|V 2 5(e)||i6||H|| L6 ||VH|| L6 + ||V^)|| L 6||VH|| 2 6 )||V 4 u|| l2 

+ (||V^)|| L 6||H|| L 6||V 2 J B|| L 6 + || 5 ( e )|| LO o||V J B|| L 6||V 2 H|| L 3)||V 4 U || L2 
+ \\g(g)h~\\B\\ L ~\\V 3 B\\ L 4\7 4 u\\ L2 (3 ' 27) 

< ||VH|| 2 2 ||V 2 H|| 2 2 + ||V 2 H|| 2 2 ||V 2 H||^- 1 +e||V 4 u||| 2 

<\\VB\\ 2 Hl \\V 2 B\\ 2 Hl +e\\V\\\ 2 L2 

< (1 + t) 5 + e||V 4 u||| 2 . 

In the same manner, it arrives at directly 

IV 7 < (l + t)- 5 + e||V 4 u|| 2 2 . (3.28) 

Substituting (j3.22j) - (!3.28l) into (13.211) . we obtain 

^ |(|V 3 0 | 2 + |V 3 u| 2 )dx + /i J |V 4 u| 2 dx< (l + t)- 5 + eo ||V 3 e||| 2 , 

which, together with (13.201) . completes the proof of the lemma. □ 


Next, we establish the inequality to recover the dissipation estimate for g. 


Lemma 3.4. Under the assumptions in Theorem \1.3l the global classical solution ( g , u , B ) of Cauchy 
problem m-m satisfies 


d_ 

dt 


j V 2 u • \7 3 gdx + Cq 


/ 


V 3 g\ 2 dx < C 7 [(l + t)“ 5 


+ ||V 3 u 



(3.29) 
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Proof. Taking V 2 operator on both hand sides of (12.11^ . multiplying by V 2 g and integrating over 
R 3 , then we have 

J(V 2 u t ■ V 3 g + \X7 3 g\ 2 )dx = J[fj,AV 2 u + (n + u)V 3 divu + V 2 S 2 ]V 3 £ dx. (3.30) 

In order to deal with the term j V 2 ut ■ V 3 gdx, we turn the time derivatives of velocity to density 
and apply the transport equation (12.1 D x . More precisely, we get 


/ 


V 2 ut ■ V 3 gdx = -jj J V 2 u • V 3 gdx — J V 2 u • V 3 gtdx 

= Tkt. f ^ 2u ' j ■ V 2 gtdx 

= -jj f V 2 ii ■ V 3 gdx — f V 2 divu • V 2 (divu + gdivu + uS7 g)dx. 


Substituting (13.31 p into (13.301) . it arrives at 
d 


dt 


V 2 u • V 3 gdx + J \\7 3 g\ 2 dx 
= J | V 2 divu| 2 dx + J V 2 divu • V 2 (gdivu + uS7g)dx + J V 2 5 2 • W 3 gdx 
+ J [/jAV 2 m + (fx + V 3 divu] V 3 gdx. 

With the help of time decay rates (13.21) . Holder and Sobolev inequalities, we obtain 
J V 2 divu • \7 2 (gdivu + uS7 g)dx 
= — V 3 divu • V(^divn + uV g)dx 


< 


< 


(||Ve||L3[|V«|| L « + ||£|M|V 2 u|| L 3 + ||n|| L6 ||V^|| L 3)||V 4 n|| L2 


l|Ve||^||V 2 n||^ + llVull^HV 2 ^ +e\\\/%\\h 


74„.ii2 


<(i+r 5 + nv 4 < 2 . 

On the other hand, just following the idea as (|3.24p ~ (|3.27l) . we have 

J V 2 S 2 ■ V 3 gdx < (1 + t)~ 5 + || V 4 n|| 2 2 + e||V 3 ^|| 2 2 

and 


J[gAV 2 u + (n + n)V 3 di .vu]V 3 gdx < ||V 4 n|| 2 2 + e||V 3 e|| 2 2 . 
Plugging (I3.33p - (I3.35D into (|3.32l) . we complete the proof of lemma. 


(3.31) 


(3.32) 


(3.33) 


(3.34) 

(3.35) 

□ 


Now, we establish optimal decay rates for the second order spatial derivatives of density and 
velocity. 

Lemma 3.5. Under the assumptions in Theorem \1.3l then the density and velocity have following 
decay rate 

\\V 2 Q(t)\\m + ||V 2 u(f)|| H i < C(1 + t)~* (3.36) 

for all t > T*(T* is a constant defined below). 
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Proof. Multiplying (13.291) by 2C ^° and adding to (13.131) . then we have 

+ C 8 J (|V 3 e| 2 + |V 3 u| 2 + \\7 4 u\ 2 )dx < C 9 ( 1 +f)" 5 , 

where £f (t) is defined as 

£%(t) = ||V 2 ,||^ + ||V 2 u||^ 1 + [ V 2 u • \7 3 gdx. 

By virtue of the smallness of £q, we have 

Cf 4 \\S7 2 ( e ,u)\\ 2 Hl < si(t) < c 10 \\v 2 ( e ,u)\\ 2 Hl . 

It follows directly from (13.371) that 

jfKt) + Y J O^ 2 + l V3 ^ 2 + l V3 “| 2 + |V 4 n| 2 )dx < C 9 (l + t)~ 5 . 

In the same manner as (12.721) . we have 

J |V 3 £| 2 dx >j^- t I \V 2 e\ 2 dx-(j^ j \\/g\ 2 dx, 

and 

\\Vuf Hl . 

Plugging (13.411) and (13.421) into (13.401) . it follows 


i^> + t 


R 

1 + t 


(|V 2 £»| 2 + | v up + |V 3 up)dx + / \X7 3 g\ z dx 


< 


< 


R 

1 +1 


{\Vg\ 2 + | V< + \V 2 u\ z )dx + (1 + t) 


\—5 


(l+f)~ 2 (l+f)"2 + (l+t) 


-5 


<(1 + i)-s. 

For some large time t > R — 1, we have 

which implies 


R 

1 +t 


< 1, 


j\V 3 g\ 2 dx < j \S7 3 g\ 2 dx. 
Combining ()3.43j) with (13.441) . it is easy to deduce 


— F 3 ( f \ , C$R ||V7 2 
df 2 ^ ' + 2(1 +t) 


V 2 (g,u)f Hl <(l + t)-*: 


which, together with the equivalent relation (13.391) . yields 

Jt S2{t) + 2 Cw(l + t) £ * {t) ~ (1 + r§ - 


(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 

(3.42) 

(3.43) 


(3.44) 


(3.45) 
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If choosing R = in (|3.45[) , it is easy to deduce 

—^ (t) ~ (1 +*) 2 > (3.46) 

for all t > T* := ^+2- — 1. Multiplying (I3.46P by (1 +7) 4 , then we have 

f[(l + t)% 3 (t)]<(l + t)~l (3.47) 

Integrating (13.4711 over [0,7], then we get 

£i(t) (1 + 1 ) 2 j 

which, together with the equivalent relation (13.391) . gives 

\\V 2 gm 2 m + \\V 2 u(mm<C(l + t)-^. 


Therefore, we complete the proof of the lemma. □ 

Finally, we establish optimal decay rate for the third order spatial derivatives of magnetic field. 

Lemma 3.6. Under the assumption of Theorem 17.51 then the magnetic field has following decay 
rate for all t >T*, 

\\V 3 B{t)\\ L2 <C(l + t)~*. (3.48) 

Proof. By virtue of (|3.9I) and applying time decay rates (|3.3I) and (|3.36l) . we obtain 

d 


-f / \V 3 B\ 2 dx + 
dt I 


\\7 4 B\ 2 dx 

< ||V(u,B)||^||V 2 (u,B)||^ 
<(l + t)-§(l + t)-i 


<(1 + t )~ 6 , 

which, together with (13.lip , gives directly 

I \W 3 B\ 2 dx +[\V 3 B\ 2 dx 


d_ 

dt J 1 1 1 + 7 

<(l + 7)- 2 ||V 2 B|| 2 2 + (l + 7)- 6 

<(l + 7)- 2 (l+7)-i+(l + 7)- 6 


(3.49) 


<(1 +t)~. 

Multiplying (13.4911 by (1 + 7) 5 and integrating the resulting inequality over [0,7], we obtain 

l|V 3 B(7)||i 2 <(l + 7)-i. 


Therefore, we complete the proof of the lemma. □ 

Proof of Theorem 11.31 With the help of Lemma 13.11 Lemma 13.21 Lemma 13.51 and Lemma 
13.61 we complete the proof of Theorem 11.31 
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3.2. Proof of Theorem \ 1-4 


In this section, we establish the time decay rates for the mixed space-time derivatives of global 
classical solutions. 

Lemma 3.7. Under the assumptions in Theorem, \1JA the global classical solution (g, u , B ) of Cauchy 
problem 42U-Q has the time decay rates 

\\^ k Qt(t)\\H 2 ~ k + l|V fc ut(i)||L2 < C(1 + t) 4 , 

\\v k B t m L * <c(i+t)-^, 


(3.50) 


where k = 0,1. 

Proof. First of all, applying the estimate (12.771) and time decay rates (11.91) . we obtain 


llVftHi, <(! + *)- 


Applying the equation (12.11) 1 . time decay rates (11.91) . Holder and Sobolev inequalities, it arrives at 

l|V 2 £>t||| 2 = || — V 2 divu — V 2 (£>divu + u • Vf?)||| 2 
< ||V 3 u||| 2 + ||V 2 (e,u)||i 3 ||V(e,u)||| 6 

12 '|v 3 mii ! 2 


+ ||(e,u)||to 

<(i +tr 


(3.51) 


Combining (13.501) - (13.511) with (12.761) . it is easy to obtain 

\\V k Qt{t)\\ 2 H 2- k <(l + t)-^, (3.52) 

where k = 0,1. Secondly, in view of the equation (I 2 .ID 9 . (j2.56[) and Holder inequality, we get 
||Vu t || 2 2 = H/uAVu + {pt + ^)V 2 divu - W 2 g + V5 2 || 2 2 

< l|V 3 u||i 2 + ||V 2 (?|| 2 2 + \\W( Ql B)\\ 2 H1 \\W 2 (u,B)\\ 2 Hl 
+ S(\\V 2 g\\ 2 L2 + \\V 2 u\\ 2 L2 + \\V 2 B\\ 2 L ,) 

<(1 + t)~ 7 u 

which, together with (|2.78|) . gives directly 

\\V k u t (t)\\ 2 L2 <C(l + t)-^, (3.53) 

where k = 0,1. Finally, it follows from (l2.1l) o. (I2.68I) - (I2.70I) . Holder and Sobolev inequalities that 
\\VB t \\ 2 L2 = ||VAH + V5 3 || 2 2 

< ||V 3 H|| 2 2 + \\V(u,B)\\ 2 h1 \\V 2 (u,B)\\ 2 l2 + ||V 2 e ||| 2 ||V J B||| 1 

< (1 + t)~ 9 2 + (1 + t)~l (1 + t)-l + (1 + t)"I( 1 + t)~ I 
<(1 

which, together with (I2.79j) . gives directly 

\\V k B t (t)\\ 2 L2 <67(1 + *)-^, (3.54) 

where k = 0,1. Combining (|3.52l) . (|3.53l) with ()3.54l) . then we complete the proof of lemma. □ 
Proof of Theorem 11.41 With the help of Lemma 13.71 we complete the proof of Theorem 11.41 
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